Sustained ferromagnetism induced by H-vacancies in graphane 
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The electronic and magnetic properties of graphane with H-vacancies were investigated using the 
quantum-chemistry methods. The hybridization of the edges is found to be absolutely crucial in 
defining the size of the bandgap, which is increased from 3.04 eV to 7.51 eV when the hybridization 
is changed from the sp^ to the sp^ type. The H- vacancy defects also influence the size of the gap 
that depends on the number of defects and their distribution between the two sides of the graphane 
plane. Further, the H-vacancy defects induced on one side of the graphane plane and placed on the 
neighboring carbon atoms are found to be the source of ferromagnetism which is distinguished by 
the high stability of the state with a large spin number in comparison to that of the singlet state 
and is expected to persist even at room temperatures. However, the ferromagnetic ordering of the 
spins is obtained to be limited by the concentration of H-vacancy defects and ordering would be 
preserved if number of defects do not exceed eight. 



INTRODUCTION 

Graphene is the carbon-based wonder material which 
has gained wide attention due to its many unique elec- 
tronic and magnetic properties. Despite the high mo- 
bility of the charge carriers in graphene resulting from 
its zero effective mass [H-Q, the absence of gap hin- 
ders its application in nanoelectronics. The magnetic 
properties of graphene arising from spin ordering of the 
localized states at the zigzag edges [J| or by the pres- 
ence of defects might facilitate its application in 
carbon-based spintronics. If the localized states occupy 
the same sublattice then they can induce the sublattice 
imbalance and according to Lieb's theorem [5] that can 
lead to the ground state being ferromagnetic [6, 7J. The 
room-temperature ferromagnetism in graphene has been 
obtained experimentally [11|. However, there are some 
issues involved in maintaining the ferromagnetic state 
whose stability depends on the concentration of the lo- 
calized states, distance between states and size of the 
graphene flakes through the size of the band gap 
Therefore, disappearance of the gap in bulk graphene 
brings some inconsistencies in applying Lieb's theorem 

i- 

Recently discovered graphane lB-l3 - hydrogenated 
graphene - has brought new impetus in the investiga- 
tion of carbon-based materials due to the predicted ad- 
vantages in its application in nanoelectronics and spin- 
tronics. Termination of the carbon atoms by hydrogens 
leads to the generation of sp^ carbon network removing 
the TT bands from its band structure thereby generating 
a gap. It was theoretically predicted that fully hydro- 
genated graphane is non-magnetic and a wide band gap 
semiconductor 15, However, H-vacancy defects in 



graphane generate localized states characterized by non- 
zero magnetic moments (each defect has ^ = 1.0/i^ po| . 
where is the Bohr magneton [Fig. j). As graphane is 
characterized by a wide gap and the value of the charge 
transfer integral {t^ ^-7.7 eV [3]) is higher than that 
in graphene {t^ ^-2.4 eV [l3|), according to the Hub- 




FIG. 1: The spin density for a single H-vacancy defect in 
graphane plotted with isovalues of ±0.001 e/A'^. 



bard model 0, 0] these should stabilize ferromagnetism 
in graphane (for example through an increase of the crit- 
ical value of the on-site repulsion term '^). 

We report here on our investigation of the electronic 
and magnetic properties of graphane and their modi- 
fication once the H-vacancy defects are introduced in 
the lattice. Our study is performed via the quantum- 
chemistry methods using the spin-polarized density func- 
tional theory with the semilocal gradient corrected func- 
tional (UB3LYP/6-31G) in the Jaguar 7.5 program [21|. 
The H-vacancies are introduced in the originally opti- 
mized structure of defect-free graphane in the chair con- 
formation (for the size of the graphane flake see Fig. ^ , 
whose carbon atoms at the edges are terminated by two 
hydrogen atoms thereby preserving the sp^ network over 
the whole lattice. 



SINGLE H-VACANCY DEFECT 

The size of the band gap of graphane nanoribbons de- 
creases exponentially with increasing nanoribbon width 
as a result of vanishing of the confinement effect [l3] sim- 
ilarly to that in graphene [isj . Therefore, for graphane 
flakes confined by the edges from all sides of size 
IOAxIBA suggesting strong confinement effect, the large 
band gap in comparison to that obtained for nanorib- 
bons U 15, l3| is expected. Defect- free graphane flakes 
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of size lOAxieA with ed ges possessing sp^ hybridization 
are found to be characterized by degenerate bands and by 
a band gap of 7.51 eV (the highest and lowest molecular 
orbitals are HOMO=-6.09 eV and LUMO=1.42 eV, re- 
spectively), as shown in Fig. n(a). For comparison we have 
examined flakes of size 18Axl6A and found a decrease 
in gap to 7.15 eV (HOMO==-5.87 eV and LUMO=1.27 
eV). The edges in the sp'^ hybridization, for which the 
edge carbon atoms are terminated by a single hydro- 
gen, possess the localized states. In this case the or- 
bital degeneracy is lifted and the gap decreases to 3.04 
eV (HOMO=-4.64 eV and LUMO=-1.58 eV). In the 



available experiment [16|, ll7|, only a transformation of 



graphene from the conductor to an insulator due to its 
hydrogenation was reported, but the size of the gap and 
the type of edge hybridization were not indicated. Since 
we found that the gap is sensitive to edge hybridization 
and can be increased from 3.04 eV to a maximum of 7.51 
eV by its transformation from sp^ to sp^ type, this is- 
sue should be the top priority for further experimental 
investigations. 

A single H- vacancy defect in the graphane lattice gen- 
erates an unsaturated dangling bond on the carbon atom 
- TT unpaired electron (perpendicular p^ orbital). More- 
over, bonding of the carbon atom carrying the defect 
with its neighbors is changed from the sp^ hybridization 
to sp^, thus providing a modification of the bond length 
from 1.55 A to 1.52 A. The perpendicular p^ orbital and 
the C-C bonds possessing sp'^ hybridization participate 
in the formation of the localized state characterized by 
an unpaired spin (see the spin density in Fig. J . There- 
fore, this localized state is spin-polarized and generates 
a defect level in the band gap (see the bands in Fig. jjb)). 
For the a-spin state, the defect level (tt) appears close 
to the valence band (HOMO=— 4.55 eV) thereby sup- 
pressing the size of the gap to A^=6.06 eV, while for the 
/3-spin state the tt* level occurs closer to the conduction 
band (LUMO=-0.84 eV) and the gap is A^=5.29 eV. 



DIFFERENT DISTRIBUTION OF THE 
H- VACANCY DEFECTS 

For several H- vacancy defects, ordering of spins of the 
localized states and the size of the bandgap are defined 
by the distance between the defects and the distribution 
of those defects between the sides of the graphane plane. 
The side dependence is related to the sublattice symme- 
try. For the chair conformation of graphane, the carbon 
atoms belonging to different sublattices are terminated 
by the hydrogen atoms from different sides of the plane. 
It was already known that in graphene [22|, when the 
localized states occupy the same sublattice and if their 
spins have antiparallel alignment, then the contribution 
of the TT states to the total energy diminishes as a result 
of the destructive interference between the spin-up and 



spin-down tails. Therefore, according to Lieb's theorem 
[5j for the localized states occupying the same sublattice 
a ferromagnetic ordering of their spins would be energet- 
ically favored, but for the states on different sublattices, 
one expects the antiferromagnetic ordering 0-0]. 

For vacancies equally distributed between both sides 
of the graphane plane (A_B-distribution), the energet- 
ically favorable spin ordering is the antiparallel align- 
ment of spins between one side (A sublattice) and the 
other (B sublattice). Thus, for even number of vacan- 
cies in the A_B-distribution, all spins are paired but the 
band degeneracy can be slightly lifted. In Fig. J^c) we 
present the bands for two H- vacancy defects in the AB- 
distribution separated by a distance d — 3aQ_Q, where 
is the length of the C-C bond. Therefore for the 



^c-C 



a- and /3-spin states, the obtained band gap of the size 
^ — 5.41 eV is defined by the energy gap between the 
TT (IIOMO=-4.56 eV) and the tt* (LUMO=-0.85 eV) 
defect levels. For odd number of defects, one localized 
state would have unpaired spin which generates an extra 
level TT for the a- and tt* for the /3-spin states. How- 
ever, when several H-vacancy defects are located on the 
same side of the graphane plane the parallel alignment of 
their spins would be preferred because they belong to the 
same sublattice (ylA-distribution) . In Fig. jjd) we show 
the energetics of the bands for graphane with two H- 
vacancy defects separated by a distance of d = 4aQ_Q in 
its triplet state. Each spin state induces a defect level in 
the valence band of the a-spin state (tt states) and in the 
conduction band of the /3-spin state (tt* states) . There- 
fore, the size of the bandgap for the a-state is defined by 
the conduction band and the defect level tt in the valence 
band, while for the /J-spin state by the valence band and 
the defect level tt* in the conduction band (Ajj=6.08 eV, 
A^=5.24 eV) that is similar to the case of the single H- 
vacancy (see Fig. J|^b)). However, the state characterized 
by antiparallel alignment of two spins (the singlet state) , 
possesses the tt and tt* defect levels for both the a- and 
/3-spin states and the size of the bandgap is defined by 
the energy gap between these defect levels, tt and tt* , i.e., 
A^ ~ A^ (for example see Fig. d[c)). 

The destructive and constructive contributions of the 
spin tails of the localized states decrease with increasing 
distance between the defects 22|. Therefore, we have 



calculated the difference in the total energy between the 
triplet and singlet states {E 



(2. 1 V 

V 2 ' 2 ^ 



depending on the dis- 
tance between the two H-vacancy defects. If E,2. 1 \ en- 
ergy is negative the ferromagnetic ordering of spins is 
energetically preferred, but an antiferromagnetic order- 
ing otherwise. We have calculated the two components: 
the energy E% l^ is considered before relaxation of the 

V 2 ' 2 ^ 

lattice induced by the presence of defects and the E,2 

^ 2 ' 2 ^ 

component after relaxation. These energies for the AA- 
and Ai?-distributions and splitting of the tt levels in the 
valence band {si — £2) are presented in Table HI For 
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FIG. 2: Energetics of the bands in graphane with sp^ hybridized edges (solid lines) and the defect levels (tt and tt*) induced 
by H-vacancies (dashed lines): (a) defect-free graphane; (b) graphane containing a single H-vacancy; (c) graphane containing 
two H-vacancies distributed between the two sides of the graphane plane (AB-distribution) and separated by a distance of 
d = 3aQ_c (antiferromagnetic spin ordering); (d) graphane containing two H-vacancies located on one side of the graphane 
plane (AA-distribution) and separated by the distance of ci = 4ac_c (ferromagnetic spin ordering). The spin density plotted 
with isovalues of ±0.001 e/A'^ is also presented. 



the ^B-distribution the distance between the two defects 
should be > aQ_Q because for d — ac-C spins of the 
two localized states are paired which leads to transfor- 
mation of a single bond in sp^ hybridization between the 
nearest-neighbor carbon atoms to a double bond in sp^ 
hybridization {a and tt bonds) of length 1.36 A. Such a 
defect forms the tt and tt* defect levels which are ener- 
getically close to the edges of the conduction and valence 
bands of graphane. Therefore, the size of the band gap 
defined by the defect levels is 6.38 eV, that is much larger 
than that for the levels formed by the non-bonded per- 
pendicular p^ orbital (A=3.71 eV in Fig. jjc)). 

A significant energy difference between the triplet and 
the singlet states was found only for d < 2aQ_Q in the 
AA-distribution for which the relaxation of the graphane 
lattice stabilizes the state with ferromagnetic spin order- 
ing, and for d < 2aQ_(~. in the Ai3-distribution. As a 
result, for the nearest location of the defects, the order- 
ing of spins is according to Lieb's theorem When the 
defects are spatially separated {d > 2a^_^) the energy 
difference between the state with a large spin number 
and the singlet state diminishes because of decoupling 
of the magnetic moments of the localized states and the 
random spin distribution with a minimum number of un- 
paired spins would be energetically favored. Therefore, 
for even number of defects the system would prefer to re- 
main in the singlet state independent of the distribution 
of defects over the sublattices, while for odd number of 
defects the triplet state is preferred. 



TABLE I: Difference in the total energy between the triplet 
and singlet states for the non-optimized plane of graphane 
with two defects E?2.i-, and after its full relaxation E.2 i\- 

(ei— £2) is the energy splitting of the tt orbitals for the relaxed 
lattice of graphane with two defects. 



distance 






(£1 - £2) (eV) 


AA-distribution (ferromagnetic ordering) 


d = 2ac-c 


-1.52 


-1.23 


2.82x10"^ 


d = 4ac_c 


-0.26 


-1.32 X 10"^ 


1.52x10"^ 


d = 6ac-c 


7.11x10"^ 


-1.15 X 10"^ 


1.53x10-2 


rf = 8 ac-c 


1.76x10"^ 


-7.93x10"^ 


6.52x10"^ 


y4_B-distribution (antiferromagnetic ordering) 


d = ac-c 


1.34 


2.98 




d = 3ac-c 


4.30x10"^ 


1.18x10"^ 




d = 5ac-c 


-1.08 X 10^" 


6.13x10"^ 





We have calculated the fluctuation of the gap size with 
increasing number of defects for the system in its singlet 
state when the size of the gap is defined by the energy 
difference between the induced defect levels, tt and tt*, 
formed by the perpendicular p^ orbitals. For the AB- 
distribution {d > «(-._(-.), the size of the gap can fluctu- 
ate from 1.2 to 3.7 eV depending on the distance between 
the defects and their locations. The change in the gap 
size is related to the degree of the broken sublattice sym- 
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metry. Moreover, with increasing concentration of the 
H- vacancies {N > 8) we found that the number of local- 
ized states can be smaller than the number of defects. 
For the ^^-distribution the size of the gap is found to 
gradually decrease from 3.75 to 2.72 eV with increasing 
number of defects from A'' = 2 to A'^ = 18. Additionally, 
with a growing number of defects a significant distortion 
of the planarity of graphane, such as buckling of the lat- 
tice inherent for the AA-distribution, was observed. 



STABLE FERROMAGNETISM 

A prediction of the ordering of the ^j4-distributed 
localized states on the graphane surface is controver- 
sial. According to results in 2J] for semihydrogenated 
graphene, i.e. graphene hydrogenated from one side 
(the so called graphone), the non- hydrogenated side 
of graphane is possessing the localized states and the 
spins all of them are ferromagnetically ordered. We be- 
lieve that this behavior is highly debatable because it 
is known for the carbon-like materials, such as diamond 
and graphitic structures, that the total magnetization 
is suppressed by increasing the vacancy concentrations, 
and particularly for graphitic structures this occurs more 
rapidly \lO]. In contrast, authors of Ref. [l^] did not con- 
sider the ferromagnetic ordering of the states in defected 
graphane at all because they found that for the vacancy 
defects located on the neighbors {d = 2aQ_Q), the spins 
are paired, indicating a nonmagnetic state. When the 



distance exceeds 2aQ_Q for which pairing of the spins can 
not occur, the interaction between the localized states 
vanish, thereby favoring the antiferromagnetic ordering 
of spins. However, we believe that pairing of the spins 
of the localized states located on the neighboring carbon 
atoms is unlikely, considering the significant distance be- 
tween the vacancies (~ 2.55 A). Moreover, pairing may 
occur only when the spins are antiferromagnetically or- 
dered (that for the states localized on the same sublattice 
is against the Lieb's theorem and it implies the for- 
mation of the bond, which again seems unlikely because 
of the large distance between the localized states (the 
typical C — C bond length in organic compounds is ^ 1.4 
A against ~ 2.55 A for the states separated by d = 2a(-,_(-, 
in graphane). Therefore, the question of ordering of spins 
of the localized states in their AA distribution is unclear 
as yet which we set out to investigate here. 

We found that for the AA-distribution of the localized 
states formed by the H-vacancies, the ferromagnetic or- 
dering of their spins is possible when the vacancies are 
placed on the neighboring carbon atoms (see £'?2 i-, for 
the ylA-distribution in Table |T|, thereby generating a 
state characterized by a large spin number (see an ex- 
ample for graphene in Ref. jl^l). Just as in [2^ it was 
noticed that increasing the distance between the vacan- 
cies leads to vanishing of tt-tt interaction between local- 



ized states resulting in the antiferromagnetic ordering of 
spins. Therefore, we investigating here the stability of the 
ferromagnetic ordering of the spins of the localized states 
placed on the neighboring carbon atoms (rf = 2aQ_Q) 
and formation of domains of defects depending on the 
concentration of the defects in domain and number of 
domains. 

We found that with increasing concentration of de- 
fects {N) the stability of the state with a large spin 
number, i.e., the difference in total energy between the 
state with a large spin number and the singlet states, 
can increase. Therefore, for two parallel lines of defects 
-3.32 eV for iV = 4 (Fig. jja) for the spin 

-5.03 eV for 



^2' 2' 



density distribution of = 4), E^^e.i-^ 
iV = 6 and E% i. = -6.56 eV for N = 8. However, for 

(. 2 ' 2 ) 

N = 8 the states characterized by the lower spin num- 
ber are close in energy to the state with the larger spin 
number {E^.e, = -0.035 eV, E%.4, = -0.065 eV and 



V 2 ' 2 ^ 



-0.12 eV) thereby destabilizing it and limit- 
ing the number of defects having ferromagnetically or- 
dered spins. A further increase of the defect concentra- 
tion (N > 8) leads to significant suppression of the energy 
difference between the ferromagnetic and singlet states. 
Thus, E'^^o.i, 



-0.25 eV for N = 10, E°^,,^ = -0.17 

eV for iV = 12 and E°^,,^ = -0.14 eV for N = 14. 

Therefore, just as for graphene 0-0], there is a critical 
value of defect concentration above which the ferromag- 
netic ordering of the spins of the localized states occupy- 
ing the same sublattice no longer exists. Another destabi- 
lizing factor for ferromagnetism in graphene containing 
many H-vacancies in the AA-distribution is the lattice 
relaxation leading to the buckling of the graphane struc- 
ture (£;.4.l^ = -1.49 eV for iV = 4 and E.s.i. = 1-29 

I 2 ' 2 / V 2 ' 2 z 

eV for TV = 8). There is also a significant decrease in 
stability when the defects are divided in groups, because 
the state characterized by antiferromagnetic ordering of 
spins between the groups is close in energy to that of the 
state with ferromagnetic ordering of all spins (see the spin 
density distribution for iV = 4 in Fig. jjb)). Finally, since 
graphane is a wide gap semiconductor and possess no 
localized states at the edges which could interact (tt — tt 
interaction) with states formed by the H-vacancy defects, 
the increasing size of the graphane flakes leading to sup- 
pression of the graphane gap should not drastically alter 
the interaction of the localized states and their ordering, 
which was the case in graphene Q. 

To summarize, for the ^^-distribution of the H- 
vacancies the size of the band gap can be tuned by 
the level of hydrogenation - the gap slowly decreases 
with increasing number of defects, while for the AB- 
distribution the size of the gap fluctuates in the range 
of 1.2 - 3.7 eV depending on the level of the broken sub- 
lattice symmetry. Moreover, formation of the H-vacancy 
defects redistributed over one side of the graphane plane 
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(a) large spin state (b) singlet state 



FIG. 3: Spin density for the ^A-distribution of H-vacancies 
(a) 8 H-vacancies (-E?8 i-, ~ —6.36 eV), (b) 4 H-vacancies in 

^ 2 • 2 J 

the singlet state (-E?4 i , = 4.81 x 10^^ eV). Spin densities are 
plotted with isovalues of ±0.005 e/A^. 

(AA-distribution) and located on the neighboring car- 
bon atoms belonging to the same sublattice will gen- 
erate a stable state characterized by a large spin num- 
ber (ferromagnetic ordering). For a better stabilization 
of this state, the reorganization of the graphane lattice 
in response to the occurrence of defects should be mini- 
mized. Deformation of the graphane lattice can be mini- 
mized for free standing graphane in the low-temperature 
regime and through interaction of graphane with a sub- 
strate under the condition that the sp^ hybridization of 
graphane lattice is preserved. If the rigidity of graphane 
on a substrate could be achieved, graphane containing 
H-vacancies on one side of the plane can even become a 
room-temperature ferromagnet, that obviously has enor- 
mous potentials for application in nanoelectronics and 
spintronics. Because each defect forms a perpendicu- 
lar orbital possessing unpaired spin and, therefore, its 
contribution to the magnetization is the magnitude 
of magnetization of such room temperature ferromagnet 
can ideally be regulated by the number of H- vacancy de- 
fects. However, the number of defects to achieve the 
stable ferromagnetism at room temperature should be 
limited because above a critical value the ferromagnetic 
ordering of spins would be unstable. 

The work was supported by the Canada Research 
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